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FOREWORD 

This  report  describes  the  work  carried  out  in  the  Solid-State 
Electronics  Laboratory,  Department  of  Electrical  Engineering  and 
Computer  Science,  The  University  of  Michigan,  Ann  Arbor,  Michigan, 
on  "Millimeter-Wave  Circuit  Analysis  and  Synthesis."  This  work 
was  sponsored  by  the  Air  Force  Systems  Command,  Avionics 
Laboratory,  Wright-Patterson  Air  Force  Base,  Ohio,  under  Contract 
No.  F33615-81-K-1U29  (Project  2002,  Task  03,  Work  Unit  92). 

The  work  reported  herein  was  performed  during  the  period 
March  1,  1981  to  April  30,  1985  by  Drs.  G.  I.  Haddad,  D.  F.  Peterson, 
R.  K.  Mains,  and  M.  E.  Elta  and  Messrs.  R.  Actis,  T.  L.  Crandle, 

J.  R.  East,  R.  K.  Froelich,  M.  Radmanesh,  P.  A.  Sanborn,  and 
D.  Yang.  The  report  was  released  by  the  authors  in  April  1985. 
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SECTION  I 


INTRODUCTION 

Hie  purpose  of  this  program  was  to  investigate  various  types 
of  microwave  and  millimeter-wave  devices  and  circuits  and  device¬ 
circuit  interactions  in  order  to  advance  the  state  of  the  art  and 
understanding  of  these  devices  and  circuits  and  their  utilization 
in  optimizing  the  performance  of  oscillators,  amplifiers,  power 
combiners  and  in  other  applications. 

Several  areas  of  investigation  were  carried  out  and  included, 
in  most  instances,  device  modeling,  device  characterization,  circuit 
design  and  evaluation,  and  device-circuit  interactions.  The  major 
areas  that  were  studied  under  this  program  and  the  important  results 
that  were  obtained  are  summarized  in  the  following  sections  of  this 
report  and  include  the  following: 

1.  IMPATT  device  design  and  optimization. 

2.  Power  combining  of  solid-state  devices. 

3.  FET  device  and  circuit  modeling  and  characterization 
including  monolithic  GaAs  integrated  circuits. 

U.  Magnetostatic-wave  device  design  and  analysis. 

Several  interim  technical  reports,  journal  articles,  and 
conference  paper  presentations,  which  describe  in  detail  various 
aspects  of  the  work  performed  under  this  program,  were  issued  and 
are  available.  These  include: 
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A.  Technical  Reports 

1.  Froelich,  R.  K. ,  "Computer  Modeling  of  Millimeter-Wave  IMPATT 
Diodes,"  Report  No.  AFWAL-TR-82-1107 ,  November  1982. 

2.  Mains,  R.  K.  and  Haddad,  G.  I.,  "Capabilities  and  Potential 
of  Millimeter-Wave  IMPATT  Devices,"  Report  No.  AFWAL-TR-82- 
llUl,  November  1982. 

3.  Yang,  D.  C. ,  "Large-Signal  Characterization  of  Nonlinear 
Two-Port  Active  Networks,"  Report  No.  AFWAL-TR-8U-1031* , 

April  198^. 

4.  Actis,  R. ,  "Lossless  Symmetric  TEM  Line  IMPATT  Diode  Power 
Combiners,"  Report  No.  AFWAL-TR-8E-1035 ,  April  198U . 

3.  Radmanesh,  M. ,  "Magnetostatic-Wave  Propagation  in  a  Finite 
YIG-Loaded  Rectangular  Waveguide,"  Report  No.  AFWAL-TR- 
8I-H7I*,  November  198U. 

6.  Peck,  D.  E.  ,  "Measurement  Methods  for  Accurate  Microwave 
Characterization  and  Modeling  of  MESFET  Chips,"  Report  No. 
AFWAL-TR-8L-1177 ,  January  1985. 

r.  Journal  Articles 

1.  Blakey,  P.  A.  and  Froelich,  R.  K. ,  "On  the  Transient  Analysis 
of  Circuits  Containing  Multiple  Diodes,"  IEEE  Trans,  on 
Microwave  Theory  and  Techniques,  vol.  MTT-31,  No.  9, 

PP.  781-783,  September  1983. 

Mains ,  R.  K. ,  Haddad,  G.  I.  and  Blakey,  P.  A.,  "Simulation 
f  OaAs  IMPATT  Diodes  Including  Energy  and  Velocity  Transport 
Koua'i'ns,"  IEEE  Trans,  on  Electron  Devices,  vol.  ED-30, 

Tit.  U  ,  pp.  1327-1338,  October  1983. 

3.  Mains,  R.  K. ,  El-Gabaly,  M.  A.,  Haddad,  G.  I.  and  Sun,  J.  P. , 
"Comparison  of  Theoretical  and  Experimental  Results  for 
Millimeter-Wave  GaAs  IMPATT’s,"  IEEE  Trans,  on  Electron 

I vices.  v  T.  ED-31,  No.  9,  pp.  1273-1279,  September  198L . 

4.  -’.abaly ,  M.  A.,  Mains,  R.  K.  and  Haddad,  G.  I.,  "Effect 
rf  Dor-in.-  Profile  Variation  on  GaAs  Hybrid  and  Double-Read 
IMPATT  Diode  P’erformance  at  60  and  9^  GHz,"  IEEE  Trans,  on 
Microwave  Theory  and  Techniques,  vol.  MTT-32,  No.  10,  pp . 

13^2-1 332 ,  Cctcfcer  1984 . 

;.  El— Gab'll  y,  M.  A.,  Mains,  R.  K.  and  Haddad,  G.  I.,  "Effects 
of  Doping  Profile  on  GaAs  Double-Drift  IMPATT  Diodes  at 
33  and  4a  GHz  ’Jsing  the  Energy -Momentum  Transport  Model," 

IEEE  Trans,  on  Microwave  Theory  and  Techniques,  vol.  MTT-32, 
No.  10,  pp.  1353-1361,  October  1981+ . 


-2- 


6.  Mains,  R.  K. ,  El-Gabaly,  M.  A.  and  Haddad,  G.  I.,  "Finite- 
Difference  Numerical  Methods  for  Solving  the  Energy- 
Momentum  Transport  Equations  in  Two-Valley  Semiconductors," 
J.  Computational  Phys.  (accepted  for  publication). 

Conference  Presentations 


1.  Yang,  D.  C.  and  Peterson,  D.  F. ,  "Large-Signal  Characteriza¬ 
tion  of  Two-Port  Nonlinear  Active  Networks,"  1982  IEEE 
MTT-S  Int .  Microwave  Symp. ,  Dallas,  TX,  June  1982.  Also 

in  Proc .  pp.  3^5-3^7. 

2.  Mains,  R.  and  Haddad,  G.  I. ,  "Simulation  of  GaAs  IMPATT 
Diodes  Including  Energy  and  Velocity  Transport  Equations," 
1983  Workshop  on  Compound  Semiconductors  for  Microwave 
Materials  and  Devices  (WOCSEMMAD),  San  Antonio,  TX, 

February  1983. 

3.  Actis,  R.  and  Peterson,  D.  F. ,  "A  Lossless  Radially  Symmetric 
TEM-Line  IMPATT-Diode  Power  Combiner,"  1983  IEEE/MTT-S  Int. 
Microwave  Symp.,  Boston,  MA,  May  31-June  3,  1983.  Also 

in  Proc .  pp.  209-211. 

As  can  be  seen  from  this  list,  significant  contributions  have 
been  made  in  device  and  circuit  design,  characterization,  and 
optimization.  The  following  sections  describe  in  more  detail  the 
results  obtained  under  this  program  in  the  various  areas. 
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4  —  9 

increased  understanding  of  the  device.  However,  the  microwave 

'■  KT  structure  poser  a  basic  problem  for  the  device  modeler.  An 
X-bumi  FKT  has  a  typical  doping  density  of  1017  cm  3 ,  with  higher 
•i")'ing  for  higher  frequency  operation.  The  characteristic  length 
it.  a  semiconductor  material  is  the  debye  length: 


Ld 


2kTe' 


a) 


where  k  =  the  Boltzmann  constant, 

T  =  the  temperature, 
c  =  the  dielectric  constant, 

=  the  electronic  charge  and 
K ,  =  the  doping. 

' ho  characteristic  time  is  the  dielectric  relaxation  time: 


T 


€ 

UqNd 


(2) 


where  y  =  the  mobility. 

The  debye  length  is  approximately  0.02  pm  and  the  dielectric 
relaxation  time  is  approximately  0.01  ps  for  an  X-band  FET.  The 
;.rollem  with  device  modeling  in  two  dimensions  is  constrained  by  the 
small  mesh  size  and  the  resulting  large  number  of  calculations.  The 
;  or  switching  problem  means  doing  a  solution  approximately  100 
: imes  per  picosecond.  Through  various  approximations  and  clever 
programming  the  number  of  calculations  car:  be  reduced  to  some  extent. 

h.2.1  Simple  Approximate  Model.  The  approach  used  in 
t  his  case  is  to  derive  a  "quasi-twc-dimer.sicnal"  solution.  The 
•oisson  ar.d  electron  transport  equations,  altered  to  account  for  the 
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SECTION  IV 


FET  DEVICE  AND  CIRCUIT  MODELING  AND  CHARACTERIZATION 

1.1  Introduction.  The  purpose  of  this  phase  of  the  program 
vif  lo  develop  appropriate  theoretical  models  for  GaAs  FET 
transistors ,  to  develop  measurement  techniques  for  experimental 
evaluation  of  the  devices  and  comparison  with  the  models,  the 
utilization  of  these  models  in  voltage-controlled  oscillator  design 
•u.d  implementation,  and  the  development  of  monolithic  integrated 
circuits  incorporating  these  devices. 

The  results  obtai  tied  on  the  various  aspects  of  this  phase  of 
the  program  are  summarized  in  the  following  sections. 

t.I  FET  Modeling.  Several  types  of  models  were  developed 
that  range  from  simplified  ones  that  lend  themselves  easily  to 
analytic  solution  to  very  extensive  two-dimensional  computer  models 
that  require  significant  computer  time  for  simulation. 

During  the  past  decade  the  FET  has  become  the  single  most 
.jrtant  microwave  solid-state  device.  It  is  used  in  small  -  and 
lar^e-signal  arnpl  i f iers ,  low-noise  amplifiers,  oscillators,  mixers, 

•! v!  a  variety  of  control  applications.  The  FET  has  been  the  suttee4 
•'  numerous  papers  leading  to  a  better  understanding  of  its  operation 
■!>;  i  potent  ini.  Early  FET  models  were  analytic  approximations  that 

1-3 

elected  the  effect  of  space  charge  in  the  channel.  With  faster 

••-r  nnt**r.-  and  better  numerical  methods,  full  solutions  of  the  two- 
;  hrer.s  ior.al  carrier  transport  and  Poisson* s  equations  became  possible. 
ul  ts  from  these  finite-difference  and  finite-element  solutions 
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N-way  nonresonant  combiners  which  are  often  associated  with  various 


bandlimiting  resistive  stabilization  techniques  for  suppressing 
non-power-producing  interactions  among  the  devices.  This  combining 
design  requires  no  such  stabilizing  scheme.  Suppression  of  undesired 
odd  mode:-:  is  accomplished  in  lossless  circuits  by  an  appropriate 
combination  of  device  and  circuit  which  provides  the  necessary 
condition  for  a  stable  combiner. 

"An  experimental  realization  of  the  simplest  example  of  a 
lossless  symmetric  combiner  design  is  presented  in  two  different 
circuit  mediums.  A  two-diode  microstrip  combiner  and  a  two-diode 
coaxial  combiner  were  developed  to  demonstrate  the  design  principles. 
The  combiners  operated  in  a  stable  nonspurious  amplifier  configuration. 
A  3-dB  bandwidth  of  l6  percent  with  h.2  dB  of  gain  was  observed  in 
the  coaxial  combiner  with  9 ^  percent  combining  efficiency." 
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SECTION  III 


POWER  COMBINING  OF  SOLID-STATE  DEVICES 

3.1  Introduction.  The  purpose  of  this  phase  of  the  program 
v:i.:  to  investigate  the  properties  of  a  new  scheme  for  circuit-level 
power  combining  of  negative-resistance  devices.  Both  theoretical  and 
experimental  work  was  carried  out  in  order  to  characterize  such  a 
scheme  and  to  determine  its  potential.  The  new  approach  utilizes  the 
properties  of  symmetric  lossless  TEM  line  combining  networks,  together 
with  the  bandlinited  negative-resistance  properties  of  two-terminal 
devices,  to  achieve  stable  combiner  designs  having  improved  band- 
widths  over  other  approaches.  IMPATT  devices  were  employed  in  an 
experimental  realization  and  evaluation  of  this  scheme. 

3.1  IMPATT  Diode  Power  Combining.  The  preceding  scheme  was 
•  cyeu  in  an  IMPATT  diode  power  combiner  and  the  results  obtained 
are  given  in  detail  in  Technical  Report  No.  AFWAL-TR-89-2035  entitled 
"Lossless  Symmetric  TEM  Line  IMPATT  Diode  Power  Combiners,"  by 

Act  is ,  A;  ril  1984.  An  abstract  from  this  report  is  included  here 
"or  corn:  .eteness. 

"T:.c  results  of  ar.  experimental  investigation  into  a  new 
nr: roach  *'  circuit-level  power  combining  of  negative-resistance 
c  vice.'  are  presented.  The  approach  uses  the  properties  of  symmet¬ 
ric  Joe cl  ess  TEM  line  combining  networks  together  with  the  bandlimited 
characteristics  of  IMPATT  diodes  to  achieve  stable  combiner  designs 
having  improved  bandwidths  over  other  approaches.  The  combining 
networks  utilized  ir.  this  investigation  fall  into  the  category  of 
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equations . 
using  these 


Results  of  large-signal  simulations  of  GaAs  IMPATls 
numerical  methods  are  presented  and  compared." 
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M.  A.  El-Gabaly,  H.  K.  Mains  and  G.  I.  Haddad,  IEEE  Transactions 


>-!i  Microwave  Theory  and  Techniques,  October  198k. 

"Doping  profile  parameters  were  varied  in  a  computer  optimization 
^.f  hybrid  and  double-Read  GaAs  IMPATT  diodes  at  60  and  9^  GHz.  The 
energy-momentum  transport  model  was  used  to  simulate  each  structure. 
..ptimum  results  for  the  various  structures  that  were  studied  are 
presented  and  compared." 

U,  "Effects  of  Doping  Profile  or.  GaAs  Double-Drift  IMPATT 
Diodes  at  33  and  UU  GHz  Using  the  Energy-Momentum  Transport  Model," 
by  M.  El-Gabaly,  R.  K.  Mains  and  G.  I.  Haddad,  IEEE  Transactions  on 
Microwave  Theory  and  Techniques,  October  1981*. 

"Experimentally  determined  doping  profiles  for  double-Read 
GaAs  IMPATT  diodes  operating  at  33  ana  1*1*  GHz  are  used  as  starting 
points  for  a  computer  optimization.  A  computer  simulation 
including  energy  and  momentum  relaxation  effects  was  used  to 
simulate  these  devices  as  the  lengths  of  the  drift  regions  and  the 
integrated  charge  in  the  doping  spikes  were  varied.  The  effects  of 
these  doping  profile  variations  on  diode  performance  are  presented." 

"Finite-Difference  Numerical  Methods  for  Solving  the  Energy- 
Momentum  Transport  Equations  in  Two-Va31ey  Semiconductors,"  by 
:\ .  K.  Maine,  M.  A.  El-Gabaly  and  G.  1.  Haddad,  Journal  of  Computational 
■ hyjlcs ,  accepted  for  publication. 

"Two  finite-difference  methods  for  solving  the  energy -momentum 
transport  equations  for  electrons  in  two-valley  semiconductors  are 
analyzed.  For  each  method,  stability  analyses  are  carried  out 
including  the  electric  field  terms  and  relaxation  terms  in  the 
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articles  are  included  here  for  completeness. 

1.  "Simulation  of  GaAs  IMPATT  Diodes  Including  Energy  and 
Velocity  Transport  Equations,"  by  R.  K.  Mains,  G.  I.  Haddad  and 

P.  A.  Blakey,  IEEE  Transactions  on  Electron  Devices,  October  1983. 

"Simulations  have  been  performed  of  GaAs  hybrid  double-drift 
IMPATT  diodes  at  t'O  and  91  GHz  using  a  transport  model  which  includes 
equations  for  the  average  per-carrier  velocity  and  energy.  These 
equations  are  obtained  from  the  second  and  third  velocity  moments  of 
the  Boltzmann  transport  equations,  respectively.  The  relaxation-time 
formulation  is  used  to  characterize  the  collision  terms.  Simulations 
were  also  carried  out  for  the  same  structures  using  the  standard  drift- 
diffusion  transport  model.  It  was  found  that  inclusion  of  the  energy- 
velocity  equations  significantly  modifies  the  predicted  carrier  transport 
behavior  and  results  in  somewhat  better  RF  performance  under  large-signal 
conditions  than  that  predicted  by  the  dri ft- di f fusion  simulation." 

2.  "Comparison  of  Theoretical  and  Experimental  Results  for 
Millimeter-Wave  GaAs  IMPATT's,"  by  R.  K.  Mains,  M.  A.  El-Gabaly, 

1.  I.  Haddad  and  <T.  P.  Gun,  IEEE  Transactions  on  Electron  Devices, 
.September  1981 . 

"Emulations  are  performed  of  GaAs  IMPATT  diodes  for  which  experi¬ 
mental  results  are  available  at  20,  33  and  11  GHz.  At  each  frequency, 
simulations  are  performed  using  both  the  drift-diffusion  approximation 
and  the  energy-momentum  transport  model.  It  is  found  that  inclusion  of 
energy  and  momentum  relaxation  effects  yields  better  agreement  with 
experiment . " 

3.  "Effect  of  Doping  Profile  Variation  on  GaAs  Hybrid  and 
Double-Read  IMPATT  Diode  Performance  at  60  and  9l  GHz,"  by 
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determine  better  design.,  for  millimeter-wave  IMPATTs  and  to  predict 
the  ultimate  potential  of  these  devices. 

Several  IMPATT  structures  were  simulated  using  both  the  drift- 
diffusion  and  energy-momentum  models.  When  the  energy-momentum 
equations  were  included  for  millimeter-wave  GaAs  IMPATT  simulation, 
primarily  three  effects  were  observed.  First,  the  induced  current 
waveforms  were  altered  from  those  calculated  using  the  drift-diffusion 
simulation  due  to  velocity  under-  and  overshoots.  Second,  an  addi¬ 
tional  delay  inthe  avalanche  process  due  to  energy  relaxation  was 
observed.  Finally,  the  precollection  process  for  precollection-mode 
devices  is  altered  by  the  inclusion  of  energy  and  momentum  relaxation 
o ; feet s • 

To  compare  the  simulation  results  with  experiment,  both  models 
were  used  to  simulate  GaAs  IMPATT  diodes  for  which  experimental 
results  are  available  at  20,  33  and  GHz.  It  was  found  that 
results  from  the  energy -momentum  simulation  agreed  better  with 
experimental  results  at  these  frequencies. 

A  series  of  studies  was  carried  out  to  find  the  effect  of 
varying  several  doping  profile  parameters  on  device  performance  at 
13,  At,  f.C  and  9t  GHz.  Since  good  agreement  had  been  obtained 
between  available  experimental  data  and  the  results  of  the  energy- 
m omentum  simulation,  this  model  was  used  in  these  studies.  The 
overall  goal  was  to  predict  optimum  doping  profiles  for  GaAs  IMPATT 
diodes  at  these  frequencies. 

The  methods  employed  in  these  simulations  and  the  results 
obtained  on  millimeter-wave  GaAs  IMPATTs  are  given  in  detail  in 
the  following  journal  articles  and  the  abstracts  from  these 


"The  model  uses  the  first  three  velocity  moments  of  the 
phase-space  transport  equation.  Terms  accounting  for  the  effects 
of  collisions  incorporate  energy  dependent  relaxation  times  and 
ionization  rates.  These  functions  of  energy  have  been  evaluated 
by  requiring  consistency  with  experimentally  determined  transport 
parameters.  Stable,  accurate  and  efficient  finite-difference 
approximations  to  the  transport  equations  were  developed. 

"The  simulation  was  used  to  identify  and  examine  the  influence 
of  various  transient  transport  phenomena,  including  velocity  overshoot 
and  undershoot,  and  the  lag  of  energy  with  respect  to  electric  field. 
The  most  important  effect  was  found  to  be  the  energy  lag.  This 
delays  carrier  injection,  thereby  improving  predicted  device 
performance. 

"The  results  indicate  that  parasitics,  rather  than  internal 
transport  mechanisms,  are  dominant  in  determining  the  experimental 
performance  attainable  from  millimeter-wave  Si  IMPATTs.  The  maximum 
frequency  of  operation  of  the  IMPATT  mode  is  subject  to  a  fundamental 
limit  associated  with  the  energy  relaxation  time.  This  limit  is 
approximately  500  GHz  for  Si  IMPATTs." 

The  new  model  was  next  applied  to  GaAs  diodes.  Here  there 
were  three  principal  objectives.  First,  it  was  desired  to  know  how 
the  new  effects  included  in  the  energy-momentum  transport  model  would 
alter  the  device  performance  from  that  predicted  by  the  drift- 
diffusion  model.  Second,  the  results  of  both  simulations  were  compared 
to  available  experimental  data  in  order  to  determine  which  could 
better  explain  experimental  results.  Finally,  doping  profile 
optimizations  were  carried  out  using  the  energy -momentum  model  to 
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modifications  was  to  include  a  more  realistic  band  structure  than 


is  typically  assumed  in  Monte  Carlo  simulations.  For  consistency, 
this  band  structure  must  be  used  when  the  integration  over  final 
states  is  performed  in  the  calculation  of  the  scattering  rates 
needed  for  the  Monte  Carlo  simulation.  Also,  particle  velocities 
and  energies  should  be  calculated  using  a  more  realistic  band 
structure  during  the  simulation. 

A  computer  program  was  developed  that  numerically  integrates 
the  scattering  rate  matrix  elements  for  all  the  important  scattering 
mechanisms  over  an  arbitrary  band  structure.  The  band  structure 
for  GaAs  obtained  using  the  pseudopotential  method  was  then  used  to 
calculate  a  complete  set  of  improved  scattering  rates  for  GaAs. 

A  new  model  was  first  employed  to  determine  the  properties 
uf  Si  IMPATT  devices  at  millimeter  wavelengths.  It  was  found  that, 
since  the  relaxation  rates  in  Si  are  very  fast,  there  was  very 
little  difference  in  the  results  obtained  from  the  two  models  up 
to  approximately  200  GHz.  It  was  therefore  concluded  that  the 
drift-diffusion  model  is  adequate  for  Si  up  to  very  high  frequencies. 
The  results  on  Si  are  described  in  detail  in  Technical  Report  No. 
AFWAL-TR-82-1107  entitled  "Computer  Modeling  of  Millimeter-Wave 
IMPATT  Diodes,"  by  R.  K.  Froelich,  November  1982.  An  abstract 
from  that  report  follows. 

"A  model  of  millimeter-wave  Si  IMPATTs  has  been  developed 
which  includes  transient  transport  effects  neglected  in  the  conven¬ 
tional  drift-diffusion  model.  The  new  model  is  based  on  principles 
of  energy  and  momentum  conservation. 


-8- 


/ 


less  numerical  diffusion  than  the  other  method.  The  problem  with 
this  method  is  that  it  becomes  unstable  in  the  limit  of  low  electric 


field  and  low  particle  energy,  so  that  it  cannot  be  used  for  device 
simulation  where  these  conditions  occur.  The  alternative  simulation 
nsec  a  modified  Lax  numerical  method  (modified  so  that  the  particle 
continuity  equations  are  conservative).  The  advantage  of  this  method 
is  that  it  remains  stable  in  the  limit  of  low  field  and  low  particle 
energy.  To  achieve  this  stability,  the  method  introduces  more 
numerical  diffusion  than  the  forward-time,  upwind-drift  method. 
However,  it  was  found  that  proper  choice  of  space  and  time  steps  can 
minimize  the  numerical  diffusion  of  the  modified  Lax  method  to 
acceptable  proportions. 

Details  of  these  investigations  were  included  in  a  journal 
article.  The  article  was  submitted  to  the  Journal  of  Computational 
Physics  and  has  been  accepted  for  publication. 

During  this  program,  the  displaced  Maxwellian  method  for  cal¬ 
culating  relaxation  times  was  investigated  extensively.  Relaxation 
times  were  calculated  for  many  different  values  of  constants  in  the 
expressions  for  the  scattering  rate  matrix  elements.  For  each  set 
of  relaxation  times,  the  transport  equations  were  solved  in  the 
uniform,  dc  limit  to  determine  whether  a  good  fit  to  the  static 
velocity-electric  field  characteristic  could  be  obtained.  For  all 
the  cases  that  were  tried,  it  was  not  possible  to  obtain  a  good  fit 
to  the  velocity-field  characteristic  at  high  electric  fields. 
Therefore,  the  displaced  Maxwellian  approach  was  abandoned. 

Instead,  it  was  decided  to  modify  the  Monte  Carlo  program 
to  be  useful  at  high  electric  fields.  One  of  the  necessary 
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i l'.  ')2';°K.  Matching  each  device  to  1-ft  circuit  resistance  gives  an 
<>:;timate  of  maximum  obtainable  pulsed  power.  Finally,  CW  and 
pulsed  performance  for  all  structures  and  materials  are  compared 
over  the  entire  frequency  range  simulated." 

A  very  comprehensive  computer  model  that  includes  energy  and 
momentum  relaxation  effects  was  developed  and  applied  to  Si  and 
GaAs  IMPATT  device  simulation  and  design.  Many  different  numerical 
techniques  for  solving  the  energy -momentum  transport  equations  were 
investigated.  These  included  both  one-step  methods,  where  the 
transport  equations  are  solved  only  once  at  each  time  advancement, 
and  two-step,  Lax-Wendrof f-type  methods,  where  the  equations  are 
solved  twice  per  time  step  advancement.  The  investigation  con¬ 
sisted  of  two  main  phases.  First,  the  numerical  methods  were 
analyzed  to  determine  the  range  over  which  the  method  should  be 
stable,  whether  the  numerical  method  accurately  models  the  actual 
equations,  and  whether  the  method  conserves  particle  number,  total 
momentum,  and  total  energy.  If  a  particular  method  was  promising, 
it,  was  implemented  in  a  version  of  the  energy -momentum  IMPATT 
simulation.  Simulations  were  then  carried  out  to  test  the  per¬ 
formance  of  the  method  in  practice. 

As  a  result  of  these  investigations,  two  different  versions 
of  the  energy-momentum  simulation  were  developed.  Both  versions 
utilize  one-step  time  advancement  schemes  since  stability  problems 
were  encountered  with  the  two-step  methods.  (Another  disadvantage 
of  the  two-step  methods  is  that  they  are  more  costly.)  One 
numerical  method  is  the  forward-time,  upwind-drift  finite-difference 
method.  This  is  the  more  accurate  of  the  two  methods  and  introduces 
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Both  models  were  employed  to  determine  the  properties  of 
millimeter-wave  IMPATTs  and,  in  particular.  Si,  GaAs  and  InP 
devices.  It  was  determined  that  the  drift-diffusion  model  is 
adequate  for  Si  devices  to  very  high  frequencies,  but  for  GaAs 
the  second  model  is  most  appropriate. 

A  comparison  of  the  results  obtained  from  these  models  with 
available  experimental  results  indicates  good  agreement.  The 
results  obtained  on  millimeter-wave  IMPATTs  are  described  next. 

2.2  Properties  and  Capabilities  of  Millimeter-Wave  IMPATTs. 
Several  technical  reports  and  journal  articles  were  prepared 
describing  the  results  obtained  in  detail.  All  of  these  publications 
are  available  and  there  is  no-  need  to  repeat  the  details  here. 

A  brief  summary  is  given,  however,  for  completeness. 

The  drift-diffusion  model  was  employed  to  determine  the 
properties  and  capabilities  of  Si,  GaAs,  and  InP  IMPATTs  at  milli¬ 
meter  wavelengths  and,  in  particular,  at  30,  UU,  60  and  9^  GHz.  The 
results  ace  given  in  detail  in  Technical  Report  No.  AFWAL-TR-82-llUl 
entitled  "Capabilities  and  Potential  of  Millimeter-Wave  IMPATT 
Devices,"  by  R.  K.  Mains  and  G.  I.  Haddad,  November  1982.  The 
abstract  from  this  report  is  included  here. 

"Theoretical  investigations  of  IMPATT  diodes  are  carried  out 
at  3C,  Uo,  60  and  9U  GHz.  GaAs,  Si  and  InP  diodes  are  simulated. 
Several  single-  and  double-drift  doping  profiles  are  considered. 
Extensive  results  as  a  function  of  RF  voltage  amplitude  and  dc 
current  density  are  presented.  Taking  thermal  resistance  into 
account,  the  expected  CW  performance  of  each  structure  is  pre¬ 
sented,  such  that  the  maximum  allowable  diode  temperature 
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SECTION  II 


IMPATT  DEVICE  DESIGN  AND  OPTIMIZATION 

Introduction.  The  objective  of  this  phase  of  the  program 
was  to  determine  the  potential  and  capabilities  of  IMPATT  devices 
for  both  pulsed  and  CW  operation  at  millimeter  wavelengths  and  to 
determine  the  effects  of  material  parameters  and  device  doping  profile 
on  the  expected  performance  of  these  devices.  In  order  to  achieve 
these  objectives,  theoretical  models  were  developed  and  utilized 
for  this  purpose.  These  represent  the  best  models  available  for 
IMPATT  device  simulations  and  can  be  divided  into  two  main  categories 
as  follows: 

1.  The  first  model  is  the  so-called  "drift-diffusion"  model, 
it  incorporates  the  static  velocity-electric  field  characteristics 
of  electrons  and  holes,  and  the  ionization  rates  are  assumed  to  be 
Instantaneous  functions  of  electric  field.  Other  than  those,  the 
model  is  fairly  general  and  gives  important  information  relative 
to  device  design. 

2.  The  second  model  is  the  so-called  "energy-momentum- 
relaxation"  model  where  intervalley  transfer  time,  energy  and 
momentum  relaxation  effects,  and  ionization  rates  as  instantaneous 
functions  of  energy  rather  than  electric  field  are  included.  This 
model,  therefore,  takes  into  account  velocity  overshoot  and  under¬ 
shoot  effects  which  become  very  important  in  GaAs ,  particularly 
at  millimeter  wavelengths. 
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two-dimensional  nature  of  the  problem,  are  solved  in  one  dimension. 
Eevoral  ana’.;,  lie  expressions  are  used  to  determine  the  geometry  of 
tl'.o  channel.  These  approximations  give  a  set  of  equations  that  can 
be  solved  rapidly.  At  the  same  time  most  of  the  two-dimensional 
device  physics  remains. 

The  quasi- two-dimensional  model  separates  a  two-dimensional 
FET  structure  into  a  gate  depletion  region  and  a  conducting  channel 
The  electron  concentration,  electric  field,  electron  velocity, 
I'-ctror.  diffusion  coefficient,  and  current  density  are  assumed  to 
depend  only  on  the  x  position  with  no  y  dependence.  This  geometry 
•tr.ii  a  volume  element  from  under  the  gate  region  are  shown  in  Fig.  ' 
I r.  this  two-dimensional  volume  element,  Poisson's  equation 

|  E  •  n  dl  =  -  -2-  (N  -  Nj  dv  (3) 

j  £  a 


(for  unit  width) 


BI+1LI+1 


t  (Np  -  Nd)Ax 


BI  +  BI+1 


where  K  =  the  average  electric  field  in  the  x-direction  on  the 
left  boundary, 

E  =  the  average  electric  field  in  the  x-direction  or.  'he 
right  boundary, 

I-i  =  tii--  average  electron  concentration  within  the  element, 
n  =  the  unit  normal  vector  and 
!Ifj  =  the  channel  doping. 

If  a  y-dire-cted  field  were  present,  it  would  not  change  the  results 


iti  Eq.  2.  The  y  competent  would  not  cross  the  integration  path  of 
Eq.  1.  In  the  limit  of  small  Ax,  Eq.  2  becomes 
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(a)  FET  GEOMETRY  AND  NOTATION.  (b)  VOLUME  ELEMENT 
GEOMETRY  AND  NOTATION. 
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The  current  continuity  equation 


0  J  •  n  dn 
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can  be  written  as 


(5) 
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where  =  the  electron  current  density  at  the  left  and  right 

contacts.  The  electron  current  density  is 


J 


v(E)N(x)  -  D(E) 


dfl(x)' 

dx 


ni  the  limit  of  small  Ax,  Eq.  8  becomes 


3(BJ)  =  3 (BN) 

3x  at 


(8) 


(9) 


Again,  any  y-directed  current  would  not  change  the  result  of  Eq.  9; 
however,  a  problem  does  exist.  The  velocity  and  diffusion  coefficient 
in  Eq.  8  depend  on  the  total  electric  field,  not  just  the  x-directed 
component.  At  low  fields  the  velocity  depends  linearly  on  the  field 
and  the  diffusion  coefficient  is  constant;  however,  at  high  fields 
the  velocity  is  constant.  The  velocity  along  the  channel  depends  on 
the  angle  of  the  electric  field.7  Equation  7  overestimates  the 
current  depending  on  the  angle  of  the  current  in  the  channel.  With 
this  overestimate,  the  amount  of  charge  in  the  channel  dipole  is 
usually  correct  within  a  few  percent  and  the  resulting  drain-source 
terminal  current  is  usually  high  by  approximately  10  percent. 
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Before  Eqs.  9  and  9  can  be  solved,  an  equation  relating  the 
voltage  along  the  channel  to  the  depletion  layer  width  is  needed. 
The  voltage  drop  across  a  depletion  layer  is  related  to  its 
width  w  by 


V 


d 


N  w2 
a 


(10) 


with  found  by  integrating  across  the  width  of  the  depletion  layer 
along  a  suitable  path.  One  possible  approximation  is  to  assume  that 
the  equipotential  lines  within  the  depletion  layer  are  parallel  to 
the  gate  as  shown  in  Fig.  2.  This  approximation  is  typically  used  in 
analytic  FET  models  and  is  correct  in  the  limit  of  long  gates.1-3 
With  this  approximation,  the  voltage  drop  across  the  depletion  layer 
depends  on  the  voltage  difference  between  the  channel  and  the  gate. 
Equation  10  can  be  used  to  find  the  depletion-layer  width  and  the 
resulting  channel  geometry  factor  B(x). 

The  problem  with  this  approximation  is  that  it  always  over¬ 
estimates  the  width  of  the  gate  depletion  layer.  This  can  be  shown 
by  referring  to  Fig.  3  (taken  from  a  two-dimensional  result).5 
Although  the  constant  voltage  lines  in  the  channel  are  approximately 
parallel,  they  bunch  together  in  the  depletion  layer  near  the  drain 
end  of  the  gate  depletion  layer.  Thus  a  sloping  straight-line 
approximation  can  be  made  for  the  equipotential  lines  within  the 
gate  depletion  layer.  The  two  geometries  are  shown  in  Fig.  2.  The 
electric  field  lines  in  the  sloping  approximation  are  semicircles. 

The  width  of  the  depletion  layer  in  this  case  is 


w'  =  w  sin  0/0  , 


(11) 
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POSSIBLE  DEPLETION-LAYER  POTENTIAL  LINES, 
(a)  PARALLEL- EQU I POTENT I AL  LINES  AND  (b) 
SLOPING-EQUIPOTENTIAL  LINES. 


- CONSTANT  VOLTAGE  CONTOUNS 

. MTf  JUNCTION  SNACE 

CHANGE  i.  ATEN  EDM 


FIG.  3  CALCULATED  VOLTAGE  DISTRIBUTION  IN  A  JFET . 


(2w  /v  =  5.0,  V  =  0,  V  ,  =  5.0  V) 
g  c  sg  sd 

(a)  CONSTANT  MOBILITY  IN  THE  SOURCE-DRAIN 


CHANNEL.  (b)  FIELD-DEPENDENT  MOBILITY  IN 
THE  SOURCE-DRAIN  CHANNEL.  (KENNEDY  AND 
O'BRIEN5) 


'.•liiTf  w  =  the  depletion  layer  width  from  Eq.  10, 

0  =  the  angle  of  the  sloping  equipotential  1 ine  • show:,  in 
Fig.  .b)  and 

w'  =  the  modified  depletion- layer  width, 
long  gates  0=0  and  Eq.  11  approaches  Eq.  10.  The  effect  <.  f 
Eq.  11  compared  to  Eq.  10  is  to  reduce  the  depletion  layer  width  for 
short  devices.  In  the  solution  of  the  equations,  both  w’  and  6  depend 
or.  the  voltage  conditions  along  the  channel.  A  Newton  iteration  is 
required  at  each  space  step  to  find  the  proper  value  of  w1, 

Equations  5  and  9,  along  with  either  Eq.  10  or  11,  can  be  used 
to  find  the  electron  and  field  distribution  along  the  channel  of  an 
••’K7  if  the  geometry,  material  parameters,  and  bias  conditions  are 
specified.  A  computer  program  to  solve  these  equations  was  written 
nnd  a  variety  of  structures  were  studied.  The  results  for  a  typical 
Si  and  GaAs  FET  are  presented  and  the  dc  results  are  used  to  describe 
small-  and  large-signal  operation. 

U.2.1a  DC  Results.  The  dc  characteristics  of  a  Si 
FET  are  shown  in  Figs.  3  and  U.  The  device  has  a  channel  doping  of 
1017  cm-3,  a  channel  height  of  0.25  pm,  a  drain-source  spacing  of 
0.5  pm.  and  a  0.2-pm  long  gate.  The  FET  quantities  are  all  per 
centimeter  of  gate  width.  Figure  4a  shows  the  drain-current  vs. 
irain-source  voltage  for  a  range  of  gate  voltages.  Under  dc  conditions 
the  channel  of  an  FET  is  charge  neutral  with  the  channel  dipole 
charge  balanced  and  with  the  fixed  doping  charge  in  the  depletion 
layer  balanced  by  the  corresponding  image  charge  on  the  Schottky- 
barrier  metal  gate.  Figure  Ub  shows  the  number  of  electrons  on  the 
rate.  Increasing  the  gate  voltage  (more  negative)  increases  the 
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Si  FET  PROPERTIES.  (a)  DRAIN  CURRENT  VS.  DRAIN-SOURCE 


VOLTAGE,  (b)  GATE  CHARGE  VS.  DRAIN-SOURCE  VOLTAGE,  AND 
(c)  CHANNEL  CHARGE  VS.  DRAIN-SOURCE  VOLTAGE. 


v  of  the  gate  rtvpletior.  layer.  The  volume  and  electron  number 

:i  I  r.o  depend  on  the  drain-source  voltage  as  shown.  The  number  of 
■  'Xi'L'bs  electrons  in  tt.e  accumulat  ion  layer  of  the  channel  charge  dipole 
is  shown  in  Kir.  ^c.  "or  low  drain-source  voltages  the  mobility  is 
a;  :  tvxim.Htely  constant  and  the  channel  dipole  is  small.  Increasing 
the  drain-source  voltage  increases  the  channel  electric  field  and 
increases  the  number  of  electrons  in  the  dipole  accumulation  region. 
Increasing  the  gate  voltage  (more  negative)  increases  the  volume  of 
the  gate  depletion  layer.  Increasing  the  gate  voltage  (more  negative) 
increases  the  size  of  the  gate  depletion  layer  which  decreases  the  size 
of  the  channel  dipole.  Thus  channel  charge  and  gate  charge  go  in 
opposite  directions  with  gate  voltage.  The  channel  charge  is 
approximately  one  tenth  the  gate  charge. 

The  dc  characteristics  of  a  GaAs  FET  with  the  same  geometry  and 
,io:ing  as  the  previous  Si  device  are  shown  in  Fig.  5.  Since  the  low 
field  mobility  in  GaAs  is  approximately  four  times  larger  than  in  Si, 
the  low  drair.-source  voltage  linear  region  resistance  is  also  four 
times  smaller.  The  GaAs  device  also  reaches  its  peak  velocity  at  a 
lower  drain-source  voltage.  Further  increasing  the  drai n-snuree 
voltage  drops  the  average  electron  velocity  and  the  current.  The 
gate  charge  is  similar  to  the  Si  results  with  slight  differences 
because  of  the  voltage  distribution  along  the  channel  and  different 
material  dielectric  constants.  The  shape  of  the  channel  charge  vs. 
drain-source  voltage  curve  depends  on  the  shape  of  the  mobility  vs. 
electric  field  curve.  In  Si  the  mobility  is  a  monotonically  decreasing 
function  of  field  so  the  charge  curves  in  Fig.  Uc  are  smoothly 
increasing  functions  of  field.  In  GaAs,  the  mobility  is  approximately 
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FIG.  5  GaAs  FET  PROPERTIES.  (a)  CHANNEL  CHARGE  VS.  DR^TN- 
SOURCE  VOLTAGE,  (b)  GATE  CHARGE  VS.  DRAIN-SOURCE 
VOLTAGE,  AND  (c)  DRAIN  CURRENT  VS.  DRAIN-SOURCE 
VOLTAGE. 
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coi.st.aiit  at  low  fields  ar.d  drops  rapidly  for  fields  above  the  peak 
velocity  field.  In  GaAs ,  the  channel  charge  is  approximately  zero 
for  peak  channel  fields  below  3500  V/cit  ar.d  then  increases  more 
raj  idly  than  the  Si  case  for  larger  fields  or  drain-source  voltages. 

The  current,  channel  charge,  and  gate  charge  information  shown 
in  Figs.  1*  and  5  were  generated  for  a  variety  of  FET  structures.  The 
information  generated  by  the  de  computer  program  is  stored  in  tables 
ar.d  forms  the  basis  of  the  small-  and  large-signal  models  discussed 
next . 

4.2.1b  Quasi-Static  Charge  Control  Model.  There  are 
several  possible  approaches  to  describing  small-  and  large-signal  FET 
operation.  The  first  approach  is  to  use  measured  data  and  nr.  equivalent 
circuit  to  model  the  performance.10  Small-signal  data  is  fitted  to 
an  equivalent  circuit  in  order  to  determine  the  desired  small-signal 
parameters  over  a  range  of  bias  conditions.  The  small-signal  equivalent 
circuit  can  then  be  used  to  predict  large-signal  performance  and  for 
circuit  design.  This  method  works  well  for  transistor  circuit 

s.-..--- -n  11 

: e s it  n . 

A  second  approach  is  to  derive  a  small-signal  equivalent  circuit 
using  an  analytic  dc  solution  as  discussed  in  References  1  through  4. 

An  extension  of  this  method  is  given  by  Mad jar  and  Rosenbaum . 1 2  The 
experimental  approach  depends  on  using  measured  data  that  limits  the 
method  to  existing  devices  and  requires  time-consuming  measurements. 

The  dc  small-signal  method  is  limited  by  the  assumptions  of  the 
underlying  dc  model. 

The  model  adopted  here  uses  information  from  the  dc  computer 
i re gram  to  define  a  small-signal  equivalent  circuit.  Since  the 
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underlying  dc  program  uses  fewer  assumptions  than  the  models  of 
References  1  through  A  and  12,  the  small-signal  model  should  also  be 
better.  Also,  most  of  the  small-signal  parameters  depend  on  measurable 
quantities,  so  the  model  can  be  used  as  a  starting  point  for  experi¬ 
mental  measurements. 

The  charge-control  model  used  is  shown  in  Fig.  6.  The  Q  and 

§ 

boxes  are  the  charge  stored  in  the  gate  region  and  in  the  channel 

region.  I  is  the  channel  current  (under  dc  conditions  the  channel 

current  is  equal  to  the  drain  current).  The  quantities  Q  ,  Q  and  I 

eg  c 

are  found  using  the  dc  program.  Under  time-varying  conditions,  the 
gate  current  is 


and  the  drain  current  is 


idM  « 


♦  wv 


A  small-signal  equivalent  circuit  can  be  derived  from  the  charge-control 

model  by  performing  a  small-signal  expansion  of  0  ,  Q  ,  and  I.  The 

eg  c 

small-signal  expression  for  Qc  is 


Q  (t) 
c 


=  Q  +  AQ  e  =  Q  + 

CO  C  CO 


3Qc  Vt  KC 
vrr  AV  e  +  —r 
3V  g  ?  •! 

g  : 


where  A0_c  -  the  small-signal  channel  charge, 

AV  =  the  small-signal  gate  voltage  and 

g 

AV^  =  the  small-signal  drain-source  voltage. 
Similarly , 
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(t)  = 


Q  +  AQ  e' 
go  g 


.  ,  3Q  .  .  3Q 

Ja)  =  <5  +  TT^  AV  eJUt  +  — ^  AV  e 

go  3V^  g  3Vd  d 


j  wt 
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j.  a 


t)  =  I  +  Al,ejuit  =  I  +-/aV  e,j“t+  ~  *V .« 

co  d  do  3V  g  3V ,  a 

g  d 


;i6) 


where  \Q  =  the  small -signal  gate  charge  and 

p 

Al  =  the  small -signal  channel  current, 
d 

Using  Eqs .  1-  through  it  in  Eqs.  12  and  13  yields 

Al  (t)  =  jwC  AV  +  jwr  AV 
f  g  g  id 


Al.(t)  =  g  AV  +  g  AV  +  jcut  AV  +  juiC  AV 
c.  in  ?  o  d  2  g  o  d 


(is; 


3Q 

£_ 

3  V 


=  the  gate  capacitance  , 


3Q 

.  = 

•t  -  3Vd  ’ 


3V 


the  transconductance  , 


3  v 


-  the  output  conductance  , 


_ c 

3V 


3Q 
_ c 

3V. 


the  output  capacitance 
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but  is  reconsidered  briefly  to  the 
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is  circuit  was  developed 
tent,  that  it  relates  to  the  calculations  present  here. 

The  operation  of  the  harmonic  combiner  relies  on  the  matching 
two  active  devices  that  are  coupled  through  a  physically  symmetric 
'.work  into  a  load.  For  this  case,  the  y  matrix  for  the  matching 
twork  and  load  is  of  the  form: 


r 

I  y  Y 

i  11  ‘12 


Y 

12 


Y 

1 1 


identical  sources  are  then  applied  at  the  inputs  of  this  two  port, 
e  eigenvectors  and  eigenvalues  of  the  system  are  given  by: 


u„ 


r 

1 — ! 

1 _ 

1 — i 

1 _ 

jkma 

II 

k  =  1  m  =  1 


k  =  2  m  =  1 


k  =  the  frequency  multiple, 

a  =  Y”/'?i 

the  number  of  identical  ports  being  combined.  It  is 
•  for  equivalent  sources,  the  voltage  on  the  load  for 

whore  m  =  1  is  restricted  to  second  harmonic  if  higher- 
rents  fire  neglected.  The  fundamental  components  from 

•  ,.f  fc-ti  v  !  y  cancel  one  another.  By  properly  designing 
■1  hi:.-  work ,  *  i.e  oscillation  of  the  circuit  can  be 
•1  :  '  ')  i.-  rr  I.-  only.  Since  for  most  cases  the  higher-order 

,  *  he  loo  :  absorbs  only  second-harmonic  power.  The 
-,g  ■  which  this  circuit  requires  is  ideally  obtained  by 
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-  .1  Ir.troduct 

i  on  . 

Tn  i  s 

T 

has< 

•  re.  ti  ■  *-  nr.-i  fabrication  of  monolithic  GaAs  integrated  cir- 
:  an:,  ::i  particular,  a  symmetric  second-harmonic  rower  combiner. 

■  circuit,  which  has  yet  to  be  implemented  mono I i th i c  al ly ,  has  a 


mendous  potential  for  extending  the  power  frequency  limitations  of 
ariety  of  active  devices.  The  goals  of  this  work  were  the  fafcri- 
ion  of  this  circuit  along  with  the  development  of  an  appropriate 
lysis  to  allow  for  the  optimization  of  circuit  performance  and  the 
rn i  o ! a :  i  \'h  of  these  techniques  tc  future  circuits  cf  this  type. 

The  results  rersented  here  describe  three  facets  of  monolithic 

■  cjrcul •  fabrication  as  well  as  some  general  characteristics  cf 

■  cl'!.;,  i  f  circuits.  A  brief  iiscussion  of  the  characteristics 

' : •; r  *u  1 1  under  consideration  is  presented  in  the  next  section. 

••  nr  foi  lowing  describe  three  areas  that  have  been  developed 
of  thi.;  work:  the  fabrication  of  GaAs  MESFETs  and  some 
■ure-d  characteristics  of  these  devices,  a  simulation  that  was 
1  ;  I  ; iet'-rr.i  no  the  harmonic  generation  efficiency  of  a  new 
tot  or  :i  i-v  struct. ur*-,  and  a  processing  schedule  that  makes  these 
levies-  structures  compatible  for  monolithic  fabrication. 

circuit  Description  and  Considerations.  The  circu 
-.r  i  which  this  segment  of  the  project,  is  directed  is  a  symmetric 
su-ncy  doubler  and  tower  combiner.  The  theoretical  basis  for 
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from  a  model  based  or.  de  and  low-frequency  parametric  measurements 


at  the  device  terminals." 

Ttolmicnl  Report  No.  AFWAI.-TR-81-1  03*t  ent.il. led  "T,urge-Signal 
' i . ; ! meter i zation  of  Nonlinear  Two-Fort  Active  Networks,"  by  D.  C. 
Yang,  April  198**.  The  abstract  of  the  report  follows. 

"The  objectives  of  this  investigation  include  the  development 
of  a  large-signal  characterization  procedure  for  nonlinear  active  tvo- 
rt  networks,  analysis  of  the  large-signal  properties  of  MESFETs , 
ar.d  establishment  of  a  large-signal  MESFET  equivalent  circuit  model. 

"An  approach  for  characterizing  the  single-frequency  properties 
of  a  nonlinear  two-port  network  is  presented.  The  characterization 
results  are  useful  in  optimizing  the  external  circuits  for  particular 
a;  :•! icationr  such  as  amplifiers  and  oscillators.  The  information 
obtained  regarding  network  properties  can  be  used  for  establishing 
Hi.-.pl:  fier  and  oscillator  design  and  stability  criteria. 

"A  GaAs  MESFET  wan  used  as  the  nonlinear  active  two-port 
network,  ar.i  its  large-signal  characteristics  were  measured.  The 
nonlinear  properties  of  the  MESFET  were  investigated  over  the  7- 
• :  i  --'Pin  money  range.  The  common  drain  configuration  was  studied 
•  o:  i  veiy  since  this  configuration  is  particularly  advantageous  for 

.  o  i  1 1  a*  -jr  1  icatioi.n  .  The  characterization  results  were  applied 
h.-sig!  ■  common  drain  MESFET  oscillators,  and  relationships 
:  ••’veer.  *  t  *••  ..utt  ut  power  and  load  impedances  were  established.  The 
_.e  j.erf  •  rmar.ee  as  a  voltage-controlled  oscillator  (VCO)  was  also 
1  :.v  igat-  i  ,  with  a  uir.rle  frequency  oscillator  built  to  verify 
urem»T.t  results. 
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the  program  is  described  in  great  detail  in  the  following  technical 
reports  and  a  brief  description  of  the  results  is  included  here. 

1.  Technical  Report  No.  AFWAI -TR-8U-H77  entitled  "Measurement 
Method  for  Accurate  Microwave  Characterization  and  Modeling  of  MESFET 
Chips,"  by  D.  E.  Peck,  January  1985.  The  abstract  from  this  report 
I'ollcws . 

"The  subject  of  accurate  measurement  and  modeling  of  three- 
terminal  chip  devices  such  as  bipolar  junction  and  field-effect 
transistors  has  been  addressed  by  others  with  varying  degrees  of 
success.  This  paper  discusses  a  two-tier  approach  for  determining 
highly  accurate  wide-band  scattering  parameters  of  chip  GaAs  MESFETs 
that  are  referenced  at  the  bonding  pads  of  the  device.  In  a  two-tier 
approach  the  first  level  of  characterization  determines  error  correc¬ 
tions  for  a  measurement  system.  The  second  level  is  characterization 
of  any  test  fixture  used  in  the  measurement  of  the  device.  The  discussion 
is  directed  at  the  second  level  of  characterization  in  which  methods 
of  oi.e-rort  untermination  and  deembedding  are  extended  to  two-port 
measurements  and  a  technique  known  as  'cold  chip’  parasitic  error  correc- 
' ior.  is  introduced  which  makes  possible  accurate  placement  of  the 
.  -parameter  reference  planes  at  the  terminals  on  the  chip  device. 

"Microwave  modeling  of  the  GaAs  MESFET  can  be  approached  from 
two  different  directions,  that  of  device  physics  and  that  of  device 
measurements.  This  paper  examines  two  different  measurement -based 
r.c.ie!  ing  methods  in  which  accurate  chip  level  S-parameter  data  is  a 
necessity.  The  first  method  uses  the  data  directly  to  determine 
element  values  of  an  equivalent  circuit  model.  In  the  second  method 
3-parameter  data  is  used  for  comparsion  to  S-parameters  predicted 
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T-V  CHARACTERISTICS  FOR  CONVENTIONALLY  DOPED  GaAs  FETs 
( v  =  's.O  x  106  cm/s  and  V,  .  =  -  l.LlTl  V) 


SOURCE 


FIG.  11  GaAs  CONVENTIONAL  FET.  (V  =  0.0  V,  V.  =  -  1.1+171  V, 

gs  bi 

DOPING  =  1.5  x  10 17  cm-3,  DEVICE  LENGTH  =0.68  um, 
DEVICE  WIDTH  =  0.11  um,  AND  GATE  LENGTH  =  0.035  urn) 


/ 


;t.ore  accurate  method  is  required.  Faricelli  et  al.13  studied 
transient  behavior  using  a  two-dimensional  transient  analysis.  Yang 
and  Peterson14  described  the  experimental  large-signal  characterization 
of  FKTs. 

A  simple  and  efficient  dc  and  quasi-static  RF  model  for  FET 
operation  was  presented.  This  model  includes  much  of  the  device 
physics  neglected  by  simpler  analytic  models,  while  at  the  same  time 
u-inr  much  lower  in  cost  than  two-dimensional  simulations.  The  model 
was  used  to  predict  the  small-signal  behavior  of  FETs.  The  model 
shows  the  importance  of  the  channel  dipole  and  the  resulting  effect 
on  the  transconductance.  Quasi-static  large-signal  operation  and 
the  resulting  limits  were  also  discussed. 

h.2.2  Computer  Modeling  of  FET  Devices.  A  two- 
dimensional  drift-diffusion  computer  model  was  developed  that  can 
he  applied  to  Si  or  GaAs  devices.  Typical  results  from  this  model 
fur  a  GaAs  FET  are  shown  in  Figs.  11  and  12.  Because  of  the  expense 
involved  Ln  using  this  model,  it  was  not  used  extensively  to 
simulate  many  devices.  It  was  used,  however,  to  compare  with  the 
approximate  model  described  in  the  previous  section  and  good  agreement 
was  obtained.  The  computer  program  was  tested  and  appears  to  be 
working  properly. 

h,3  Experimental  Characteri zation  of  GaAs  FETs.  Techniques 
were  developed  for  experimental  characterization  of  GaAs  MESFETs 
find  for  measurement  of  small-  and  large-signal  equivalent  circuit 
parameters.  This  is  extremely  important  for  determining  the  properties 
of  these  devices  and  their  utilization  in  amplifiers,  oscillators, 
and  in  other  applications.  The  work  performed  under  this  phase  of 

-39- 


/ 


i 


LO  LARGE-SIGNAL  QUASI-STATIC  OPERATION  OF  A  GaAs  FET 


time.  The  drain  current  vs.  time  is  shifted  with  respect  to  the  gate 
voltage  by  the  t  delay  discussed  previously.  The  delay  time  for  this 
case  is  approximately  5  ps.  The  current  is  slightly  clipped  at  higher 
gate  reverse  bias  because  of  the  transconductance  reduction  there. 

The  channel  charge  vs.  time  is  shown  in  Fig.  9d.  Two  effects  acting 
together  produce  the  double  peak  curve  shown.  The  first  peak  corres¬ 
ponds  to  the  minimum  gate  reverse  bias  and  maximum  channel  dipole 
volume.  The  second  peak  occurs  when  the  drain-source  voltage  is 
approaching  a  maximum,  producing  a  larger  electron  density  in  the 
dipole.  The  drain-source  voltage  and  operating  load  line  are  shown 
in  Figs.  9e  and  9f. 

A  set  of  curves  for  the  GaAs  device  of  Fig.  3  is  shown  in  Fig. 

10.  For  this  example,  R^  =  1  Q,  =  0.7  and  w  =  2  V.  The  gate 

is  biased  with  a  -  1  V  gate  voltage  and  a  1  V  RF  voltage.  The  three 

sets  of  curves  correspond  to  frequencies  of  20,  50  and  100  GHz.  The 

combination  of  R  and  the  gate  capacitance  lowers  the  RF  voltage 

g 

across  the  gate  terminal  with  increasing  frequency.  This,  in  turn, 
reduces  the  other  quantities  shown  in  Fig.  10. 

The  performance  of  the  GaAs  device  is  similar  to  that  of  the 
Si  device,  with  small  differences  because  of  the  shape  of  Q^.  This 
similarity  points  out  the  basic  limitation  of  quasi-static  models, 
or  any  model  that  derives  large-signal  information  on  the  basis  of 
small-signal  calculations  or  measurements.  The  quasi-static  model 
assumes  that  conditions  in  the  channel  are  described  by  the  stored 
charge  and  the  channel  current.  However,  depending  on  the  bias  point, 
RF  voltage,  and  RF  frequency,  the  channel  conditions  can  be  altered 
by  the  current  flowing  through  the  gate.  Under  these  conditions  a 


o:.ar.»V8  the  .amount  of  charge  in  the  gate  depletion  layer,  resulting 
i:  nr.  additional  gate  current  Jwt  AV  . 

The  small-signal  FET  parameters  depend  on  the  dc  operating 
;  chosen.  Typical  small-signal  parameters  for  the  devices  depend 

. :  t ho  parameters  ir.  Figs.  U  and  5-  This  model  is  a  model  for  the 
intrinsic  KET.  Additional  elements  would  be  required  to  describe 
parasitic  losses  and  the  effects  of  bondvires. 

Equations  12  and  13  can  also  be  used  to  describe  large-signal 
p.asi-stat ic  operation.  Expanding  the  time  derivatives  in  these  two 
•  ■  -j  .at  i  r.;  .-ri  ves 


3Q  3V  3Q  8V, 
1  (t)  =  — &  — &  +  1 

«  3Vg  3t  3Vd  3t 


(20) 


ar/: 


9Q  3V  8Q  av 

I"tl  ■  K 


(21) 


These  equations,  along  with  information  on  the  circuits  connected  to 
the  drain  and  gate,  describe  the  large-signal  operation  under  a 
variety  cf  drive  conditions. 

A  very  simple  bias  circuit  is  shown  in  Fig.  8.  The  FET  is 

biased  with  1-fi  gate  and  drain  resistance,  a  drain  supply  voltage 

of  3  V  and  a  dc  gate  bias  of  -  x . o  V.  Figure  9  shows  the  operation 

of  a  Si  FET  shown  in  Fig.  4  for  a  1.5-V  RF  gate  voltage  at  25  GHz. 

The  gate  voltage  is  shown  in  Fig.  9a.  The  combination  of  R  and  the 

8 

gate  capacitance  lowers  the  applied  peak-to-peak  voltage.  The  phase 
is  also  delayed.  The  corresponding  number  of  gate  electrons  vs. 


time  is  shown  in  Fig.  9b.  Figure  9c  shows  the  drain  current  vs. 


:.M'  sir.all-r.iri.ru  •  luivn’.  er.t  circuit  is  shown  in  Fig.  7-  The  sisal  1- 
.  -I  i  quo:  ,r. .1  ’  i  from  the  slopes  of  channel  charge, 

sat  e  charge,  a:.;  •  .rr--:.t  curves  giver,  in  Figs.  '•  and 

'  ar.a  r  are  *y:  .  ."il  smal  1-sirr.al  quantities.  Most  analytic 

c. .  U‘.c  ...c*  . :  era  _  negiect  eras r.e  1  cnar.ro.  Howev’-r,  tr.e  :  cc  _u." 

of  the  char.nei  choivv  ny  the  drain-source  voltage  lead;'  to  output 
capacitance. 

The  rorsconauctarce  term  in  the  equivalent  circuit  is  modified 

iy  ,’vt  .  t  corresponds  to  the  charging  and  discharging  of  the  chanr.e 

charge  dipole  as  the  channel  region  volume  is  pumped  by  the  gate 

depletion  layer  and  gate  voltage  [x  =  (3Q  /3V  )].  This  charging  and 

2  c  g 

discharging  produces  an  additional  channel  current  ju)i2  lagging  the 
normal  current  that  is  due  to  the  combination  of  the  transconductance 
and  the  gate  voltage  [lagging  since  x?  is  a  (-)  quantity].  A  typical 
•equivalent  circuit  has  a  gate  transit-time  term  e  ')UT  T  to  account  for 
the  measured  phase  delay  of  the  drain  current.  For  a  small  transit 
t  irr.e  the  transconductance  term  becomes 

goe  "  JwTT  ~  g0{l  "  ,3wtT^  ‘ 

The  : resent  equivalent  circuit  shows  that  the  measured  phase  delay  is 

also  due  to  channel  dipole  conditions,  x  depends  on  the  amount  of 

charge  in  the  channel  dipole  for  various  bias  conditions.  However, 

the-  characteristics  of  the  charge  dipole  are  approximately  constant 

iver  a  range  of  gate  lengths,  sc  x  should  not  depend  on  gate  length. 

2 

The  feedback  term  x  describes  the  coupling  between  the  drain - 
source  voltage  and  the  gate  charge .  Changing  the  drain-source  voltage 


-32- 


i 


fabricating  the  entire  circuit  monolith ically  in  a  physically  symmetric 
structure.  The  fabrication  techniques  that  are  describee  in  the 
following  sections  have  all  been  designed  with  this  factor  iri  mir.d. 

1.1.3  Results.  In  the  course  of  this  investigation 
several  tasks  were  performed  and  are  discussed  in  the  following 

] .  OaAs  MESFET  fabrication.  A  variety  of  processing  techniques 
for  the  fabrication  of  GaAs  MESFETs  were  explored.  They  were  used  to 
t  reduce  structures  that  were  evaluated  for  ac  and  extrapolated  for 
microwave  performance.  These  show  good  performance  and  the  results 
are  presented  here. 

1.  Varactor  diode  simulation  and  fabrication.  Efficient 
harmonic  genei  .cion  requires  a  high  degree  of  nonlinearity.  A 
nearly  lossless  method  for  providing  such  a  nonlinearity  is  provided 
by  the  varactor  diode.  Typically,  the  desired  structure  for  these 
devices  is  a  high -low-high  doping  profile.  A  new  device  structure 
that  utilizes  ion  implantation  to  produce  a  vertical  profile  of  this 
sort  was  simulated  and  snowed  excellent  characteristics. 

3.  Monolithic  fabrication.  The  fabrication  of  uevices,  such*  as 
MESrKTs  and  varactors,  monol ithically  requires  a  means  of  providing 
a i fferir.r  electrical  characteristics  on  a  single  substrate.  The  method 
fv  r  implementing  this  is.  selective  implantation  into  a  semi -insulating 
substrate.  This  provides  the  required  electrical  characteristics 
with  ■ i  minimum  of  process  technicality  and  a  great  deal  of  design 
:‘l  r,x  ibil  1  ty  . 

Mach  of  these  tasks-  is  discussed  in  more  detail  in  the  following 


i 

1+.4.U  MESr'ET  Fabrication.  Table  1  gives  a  skeletal 
outline  of  the  process  used  to  fabricate  the  GaAs  MESFETs  that  are 
discussed  here  as  well  as  some  characteristics  of  each  processing 
step.  Figures  13  through  15  show  measured  data  from  FETs  fabricated 
on  ion-implanted  low-noise  material.  When  the  methodology  described 
by  Fukui16  is  used,  the  maximum  available  power  can  be  determined  in 
a  straightforward  manner  from  these  plots.  The  gate  capacitance  vs. 
applied  voltage  curve  given  in  Fig.  15  and  the  transconductance  that 
can  be  determined  readily  from  Fig.  13  make  possible  a  simple  calcula¬ 
tion  of  the  unity  gain  frequency.  The  FETs  fabricated  on  this  material 
showed  good  characteristics  and  typical  unity  gain  frequencies  of 
approximately  8  GHz.  Maximum  available  power  for  typical  .devices  on 
>  this  material  was  also  calculated  and  typical  values  of  approximately 

30  mW  for  100-jjm  gate  structures  were  obtained.  A  significant 
improvement  in  both  power  and  frequency  performance  of  these  devices 
is  expected  for  substrates  with  thicker  more  heavily  doped  active 
j  layers.  The  final  design  of  the  doubler  circuit  will  necessarily 

include  the  consideration  of  the  characteristics  of  the  active  channel 
;  and  circuit  considerations  that  limit  circuit  power  consumption. 

^.4, 5  Varactor  Modeling  and  Process  Sequence.  A  nearly 
j  ideal  nvthc-i  for  generating  harmonic  power  is  provided  by  the 

varactor  diode.  The  application  of  a  sinusoidal  voltage  to  this 
ic vice  rives  a  time  variant  capacitance  that  results  in  deviations 
from  the  purely  sinusoidal  current  that  would  be  obtained  across  a 
constant  capacitance.  To  produce  a  large  harmonic  component  it  is 
genera] ly  desirable  to  maximize  the  variations  in  capacitance  vs. 
time.  The  general  characteristics  that  make  for  efficient  harmonic 
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TYPICAL  MEASURED  DRAIN  I-V  CHARACTERISTICS  FOR 


LOW-NOISE  FET. 


(volts) 


V 


9S 


)F  FORWARD-BIASED  SCHOTTKY  DIODE 


ARASITIC  RESISTANCES  (INCLUDES 


STANCE). 


0  .5  1  1.5  2  2.5  3  3.5  4  4.5  5 

VOLTRGE  (volts) 


.  15  CAPACITANCE-VOLTAGE  PROFILE  FOR  TYPICAL  SCHOTTKY  DIODE 


ON  LOW-NOISE  MATERIAL. 


generation  can  be  pointed  out  by  consideration  of  an  unrealistic 


r 


> 


device  that,  nonetheless  brings  out  the  salient  features.  This 
structure  consists  of  an  impulsive  doping  spike  at  the  surface  to 
which  the  Schottky  contact  is  made.  This  causes  the  zero  bias  deple¬ 
tion  width  to  approach  zero  and  allow  for  a  subsequently  large  variation 
in  the  capacitance  that  is  inversely  proportional  to  the  depletion 
width.  Beyond  the  zero  bias  depletion  width  it  is  optimal  to  have 
intrinsic  material  so  that  the  variation  in  depletion  width  with 
respect  to  voltage  is  maximized.  The  depletion  width  as  a  function 
of  voltage  for  this  case  is  given  by: 

U/2 


2eV 


qNT 


(23) 


where  V  is  the  reverse  voltage  and  is  the  doping  density.  Conse¬ 
quently,  the  variation  of  w  with  voltage  is  given  by: 


3w 

3V 


r  ni/2 
£ 


2qNDV 


(21) 


It  is  readily  apparent  from  this  expression  that  the  maximum  variation 
in  capacitance  occurs  at  small  reverse  voltages  and  low  dopings. 

In  practice,  this  type  of  structure  is  limited  by  a  variety  of 
considerations.  These  include: 

1.  The  thickness  of  the  barrier  through  which  electrons  must 
tunnel  to  contribute  to  current  flow  is  controlled  by  the  doping  at 
t ! l- -  surface.  In  order  to  avoid  a  deleterious  tunneling  current  it 
is  necessary  to  reduce  the  doping  at  the  surface. 

2.  The  peak  field  in  a  varactor  diode  occurs  at  the  interface 
between  the  metal  and  the  semiconductor.  If  it  is  assumed  that  the 
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impulse  cl'  doping  is  replaced  by  a  highly  doped  region,  the  field  in 
the  highly  doped  region  is  given  as: 


E(x) 


E 


1 


max 


where 


w  =  the  zero  bias  depletion  width  =  /2eV  /qN  , 
o  B  D 

V  =  the  barrier  height  and 


E  =  2 V  /w 
max  B  o 

at  zero  bias.  If,  as  suggested  by  the  discussion  of  the  idealized 
structure,  the  zero  bias  depletion  width  extends  across  the  highly 
coped  region,  the  width  and  doping  of  this  region  will  be  the  primary 
determinant  of  the  breakdown  voltage  of  the  structure.  It  is 
desirable  for  power  and  tuning  considerations  that  the  breakdown 


voltage  be  maximized.  It  is,  therefore,  desirable  to  maximize  the 
breakdown  voltage  arid  thus  minimize  the  doping  at  the  surface  or  at 
least  extend  the  length  of  the  highly  doped  region. 

i.  Parasitic  resistance  through  the  lcwly  doped  region  detracts 
•’rum  tno  overall  efficiency  of  the  device.  This  must  be  reduced  by 
‘ he  introduction  of  a  heavily  doped  buried  layer  as  well  as  by 
inerertr.ing  the  doping  of  the  intrinsic  layer.  It  is  important  tc 
note  that  the  biggest  variation  in  capacitance  occurs  as  the  depletion 
region  varies  around  tne  zero  bias  condition.  It  is  therefore  more 
: .oli  ful  to  have  a  doping  profile  that  increases  toward  the  buried 
highly  doped  layer. 

These  considerations  dictate  that  a  criterion  for  evaluating 
varactors  as  harmonic  generation  sources  must  relate  the  relative 
importance  of  harmonic  generating  efficiency  and  resistive  loss. 
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.'tri;  gi.1  ibsvard  ana .  ysis  that  results  from  the  assumption  that  the 
varactor  is  »  xci  by  a  sinusoidal  voltage  source  gives  a  figure  of 
•  hi-  varact  ■■■•r.  for  a  given  load,  the  loss  of  the  varactor 


2 


I  !  1  „  are-  rms  values  cf  the  total  and  first-harmonic  current, 

•  Iv--;  v  ,  -mu  r:  is  the  series  resistance  of  the  varactor.  This 

■  .  -  ,  :  it  minimized  for  peak  performance  and  can  be  used  as 

r  :■  ;•  c  mparisoi.  between  two  device  structures.  The  general 
tics  ie scribed  previously  are  readily  obtainable  by  a 
:■  .  •  r  the  vertical  doping  profile  of  a  properly  designed 

.  :  n;  lai.ta4  I  n.  An  ion  implant  simulation  based  on  the  Pierson  IV 

: ; .: :  r  ; ;  i .  was  implemented.  This  was  used  to  evaluate  the  capacitance 
vs.  v  ltng"  Jit.  J ,  consequently,  the  harmonic  generating  capability  of 
•t  number  of  device  structures.  Fourier  components  of  the  current 
wave  form  were  evaluated.  In  addition,  the  vertical  resistivity  through 
-he  lowly  doped  region  and  the  sheet  resistivity  of  the  buried  layer 
w-u'c  evaluate-:  sc  that  the  series  resistance  of  a  particular  layout 
configuration  could  be  determined.  This  makes  a  comparison  of  two 
jeviee  structures  possible  using  the  method  described.  Figures  16 
and  17  show  a  typical  implant  profile  and  the  associated  capacitance¬ 
's  it.  age  characteristic.  Figures  id  through  20  show  an  RF  voltage 
the  response  of  this  device  structure  to  this  voltage. 

Tut-  Fourier  coefficients  displayed  on  the  plot  show  a  good  harmonic 
generating  capability  while  maintaining  a  reasonably  lew  parasitic 
resistance.  The  fabrication  procedure  for  these  devices  is  similar 
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■HI.  16  IMPLANT  PROFILE  FOR  ENHANCED  HARMONIC  GENERATION 
(Si  INTO  GaAs ) . 
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CRPACITRNCE  (X1Q 


o 


VOLTAGE 


G.  17  CAPACITANCE-VOLTAGE  CHARACTERISTIC  FOR  THE  PROFILE  OF 
FIG .  16.  RV  AND  RS  CAN  BE  USED  TO  CALCULATE  THE 
PARASITIC  RESISTANCE  FOR  A  GIVEN  LAYOUT  CONFIGURATION. 
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CRPRCITRNCE  (X 


CURRENT 

-1400. OQ  _  600.00  2600.00 


FUND.=  9.276E+02 


I r, .  )  CURRENT  PRODUCED  BY  VOLTAGE  AND  CAPACITANCE  WAVEFORMS  OF 

FIGS.  18  AND  19.  THE  BARELY  VISIBLE  BUMP  ON  THE  AXIS  AT 
’.-5  ns  IS  THE  NEGLIGIBLE  CONTRIBUTION  OF  THE  FORWARD- 
BIASED  "DIODE"  CURRENT. 
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■  '  ;  hat  of  ’  .  ciscussed  previously  with  the  except! or.  of  an 

"'on  hunt  allows  tj-.tr  hti.ic  contacts  to  he  formed  on  the  buried  layer. 
o.<-  fu erica'  io::  :  equence  is  outlined  in  Table  1. 

■■i.l.r  Monolithic  Fabrication  Sequence.  The  device 
structures  Jirieussed  previously  can  be  made  monolithically  compatible 
by  selectively  implanting  a  GaAs  wafer  with  the  desired  profiles  on 
d 1 -'crept  regions  of  the  substrate.  This  can  be  done  by  masking  the 
vu'er  with  a  dielectric  layer  and  patterning  this  dielectric  to 
er  -.-jt.e  tubs  of  FET  and  varactor  material  on  the  substrate;  by  then 
using  a  mesa  etch  to  isolate  the  devices,  the  lithographic  requirements 

■  a’  tin:  implantation  are  minimized. 

It'  a  methodology  similar  to  that  used  to  predict  the  doping 
; ’mobile  in  the  substrate  is  used,  the  penetration  of  the  dopant  into 
t::c  masking  material  can  be  calculated.  In  this  case,  however,  the 
quantity  if  interest  is  the  number  of  accelerated  ions  that  will 
: -er.rt  rut  e  t  n rough  a  giver,  thickness  of  the  dielectric  mask.  Figure 
shows  a  plot  of  this  residual  fluence  for  the  varactor  profile 
: recent-  :  previously,  where  the  material  being  considered  is  a  silicon 
nitride  murk  layer.  The  effectiveness  of  a  mask  layer  can  be  deter- 
mL  .-•wii  ly  by  r.ue'n  a  plot.  As  an  added  safeguard,  the  isolation 
■v  •'  ‘hut  i performed  after  the  implants  have  been  performed  will 
-r.i-v-  any  dtutisticnlly  insignificant  number  of  dopant  atoms  that 
uu Yu  mat:  age  a  to  penetrate  into  tile  GaAs  substrate.  This  process 
so  picr.ce  all  .wc  for  ease  and  process  flexibility  that  is  unavailable 
:’r  it  ■-*  o  r  techniques.  A  brief  outline  of  the  complete  fabrication 


is  given  it.  Table  2. 
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Table  1 


FET  and  Varactor  Fabrication  Steps 


Process  Step 


Varactor 


Orientation  etch  Deep  etch  performed  on 
small  portion  of  wafer 
to  determine  crystal 
orientation. 


Same  as 


Ohmic  contact 


Not  reauired. 


Ohmic  contact 
format  ion 


Mesa  etch 


Schottky 
contact 
format,  ion 


Ni-Ge-Au:  300  X.-U00  K- 
2800  %.  Alloy:  2  min  at 
h90°C.  Typical  R  =1.5 
x  10“ 6  a/cm2.  C 

NH40H:H20:H202-  20:973:7 
Etch  rate  ~  33  X/s. 
Probe-surface  breakdown 
used  to  determine  pene¬ 
tration  through  active 
layer . 

Ti-Au:  700  t  -  2000  A. 
Lithography  schemes: 

AZ  resist  =  -0.75  pm 
resolution,  poor  lift¬ 
off,  AZ  resist  with 
chlorobenzene  =1.5  pm 
resolution,  good  lift¬ 
off,  AZ  resist  =  w/Si3N4 
'0.75  pm  resolution, 
good  liftoff. 


Etch  through  surface 

n-layer  so  contact  can  be 
made  to  buried  n+  layer. 
In  situ  determination  of 
dopant  concentration 
determined  by  probe- 
surface  breakdown. 

Same  as  FET. 


Same  as  FET. 


Same  as  FET. 


1 


RESIDURL  FLUENCE  IX10-9) 

oQ.OQ  2200.09  4400.0  6600.26  8800.34 


FIG.  21  RESIDUAL  FLUENCE  VS.  DISTANCE  INTO  A  Si.N,  SUBSTRATE 

3  4 

FOR  THE  ENERGIES  AND  PLUENCES  GIVEN  IN  FIG.  lb. 
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Table  2 


Monolithic  Fabrication  Sequence 


Process  Step 

Orientation  etch 
and  alignment 
marking 


_ _ Process  Description _ 

Deep  etch  performed  on  GaAs  wafer  to  determine 
crystalline  orientation  and  provide  initial 
alignment  reference. 


Implant  FET  profile  Implant  FET  regions  of  wafer  using  Si  N  mask. 


Implant  varactor 
profile  and  anneal 


Implant  varactor  regions  of  wafer  masking 


with  SijN^ 


Varactor  ohmic  etch  Etch  varactor  ohmic  pattern  so  contact  can  be 

made  to  n+  buried  layer. 

ohmic  metallization  Ni-Ge-Au  as  described  in  Table  1  onto  etched 

region  for  varactors  and  surface  for  FETs. 


Isolation  etch 


Schottky  contact 
formation 

Lumped-element 
metal 1 ization 


As  described  in  Table  1,  but  deep  enough  to 
isolate  both  FET  and  varactor  profiles. 

Forms  Schottky  contact  for  both  varactors  and 
MESFETs  as  described  in  Table  1. 

Forms  interdigitated  capacitors,  loop  inductors, 
or  high  impedance  line  sections;  could  be 
included  with  Schottky  metallization  for  a  well- 
developed  process. 


/ 


U.4.7  Conclusions.  Monolithic  fabrication  of  a  harmonic 


power  combiner  was  explored  and  appears  to  be  a  feasible  method  of 
extendi  nr  the  inherent  power  frequency  limitations  of  a  variety  of 
solid-state  active  devices.  GaAs  MESFETs  were  fabricated  and  tested 
for  use  in  this  circuit.  Varactor  diode  structures  utilizing  an  ion 
implantation  profile  to  produce  enhanced  harmonic  generation  capability 
w»re  modeled  and  show  promise  as  harmonic  sources.  Finally,  a  process 
sequence  that  makes  these  devices  monolithically  compatible  was 
presented. 

With  these  results  in  mind,  monolithic  fabrication  of  the  doubler 
circuit  should  be  feasible.  Implementation  of  the  circuit  will  require 
some  developmental  work  on  the  most  efficient  configuration  for  the 
MtoFETs  as  well  as  the  evaluation  of  feedback  capacitors  and  possibly 
loot  inductors.  These  can  probably  be  fabricated  in  a  single-level 
metal lization  that  has  been  provided  for  in  the  process  sequence  pre¬ 
sented  earlier.  Once  this  circuit  has  been  fabricated  monolithically, 
the  principles  presented  here  should  make  the  fabrication  of  higher- 
order  combiner  circuits  feasible.  This  should  allow  for  greatly 
improve d  performance  from  a  variety  of  device  structures. 
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MAGNETOSTATIC-WAVE  DEVICE  DESIGN  AND  ANALYSIS 

The  purpose  of  this  phase  of  the  program  was  to  carry  out  a 
theoretical  investigation  concerning  the  magnetostatic  modes  of  wave 
propagat i on  in  a  magnetized  YIG  slab  in  a  waveguide.  An  analysis 
was  developed  for  a  finite  YIG  slab  and  the  properties  of  magneto¬ 
static  modes  of  wave  propagation  were  obtained  as  a  function  of 
various  parameters  including  the  slab  dimensions  and  its  location 
in  the  waveguide  structure  as  well  as  magnetic  field. 

The  results  obtained  from  the  study  were  described  in 
Technical  Report  No.  AFWAL-TR-8U-II7U  entitled  "Magnetostatic-Wave 
Propagation  in  a  Finite  YIG-Loaded  Rectangular  Waveguide,"  by  M. 
R.'idmunesh ,  November  1984.  The  abstract  of  this  report  follows. 

"71  e  objective  of  this  investigation  is  to  study  magnetostatic- 
wave  propagation  in  a  rectangular  waveguide  loaded  with  a  finite  YIG 
..•lab,  to  understand  the  effect  of  physical  parameters  on  the  disper¬ 
sion  uni  'reut:  time-delay  characteristics,  and  to  ol  tain  better  design 
•.•riu-ri-'!  for  magnet  os  tn  tic-wave  devices. 

"An  analysis:  was:  carried  out  for  two  ba.  .c  geometrical  c<_n- 
rut  1  ns:  (.)  a  YIG  slab  of  arbitrary  thickness  placed  inside  a 
waveguide  in  contact  with  each  sidewall  and  (2)  a  thin  YIG  slab  of 
'"ini  to  width  placed  symmetrically  inside  a  waveguide  suet;  that  there 
is  an  ejunl  air  gap  between  the  slab  and  each  sidewall.  i- or  each  of 
these  basic  geometries,  the  magnetic  bias  field  was  assumed  to  be 
either  parallel  or  normal  to  the  slab  in  the  transverse  plane  of 
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used  to 


the  waveguide.  For  Case  1 ,  the  mode  analysis  technique  was 
obtain  the  dispersion  relations  while  in  Case  2  the  integral  equation 
method  was  used  to  derive  the  approximate  dispersion  relations. 

"Numerical  simulations  were  carried  out  for  various  geometrical 
confi purations.  Sample  numerical  calculations  of  the  propagation  con¬ 
stant,  and  delay  characteristics  are  presented  in  the  frequency  range 
of  U  to  10  GHz..  These  numerical  results  were  obtained  as  a  function 
of  the  geometry  and  magnetic  field.  The  results  obtained  from  the 
numerical  simulation  will  be  important  in  the  design  of  magnetostatic- 
wave  devices. 

"The  general  formulation  presented  in  this  study  contains  all 
the  information  given  by  the  degenerate  cases  previously  published. 

The  special  cases  of  interest  are  obtained  simply  by  permitting 
some  of  the  dimensional  parameters  to  take  on  their  limiting  values." 
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SECTION  VI 


CONCLUSIONS 

Several  important  contributions  were  made  during  the  course 
of  this  program.  These  were  described  in  the  body  of  the  report  as 
well  as  in  detailed  technical  reports,  journal  articles  and  conference 
presentations.  Significant  contributions  were  made  in  the  following 
areas : 

1.  Very  comprehensive  and  excellent  computer  models  were 
developed  for  IMPATT  device  simulations.  These  models  were  employed 
to  compare  with  experimental  results,  where  the  agreement  was  found 
to  be  very  good,  and  for  optimum  device  design  and  prediction  of  the 
potential  and  capabilities  of  these  devices  for  millimeter -wave  power 
generation. 

2.  Both  simple  and  approximate  as  well  as  comprehensive  computer 
models  were  developed  for  GaAs  MESFET  simulation.  These  models  were 
employed  to  determine  small-  and  large-signal  equivalent  circuits 

cor  these  devices  and  their  dependence  on  material  parameters  and 
device  geometry. 

3.  New  measurement  techniques  were  developed  for  characterizing 
i-’HT  devices  and  these  techniques  were  employed  to  measure  the  small- 

large-signal  equivalent  circuit  parameters. 

4.  New  ideas  were  developed  and  analyzed  for  power  combining  of 
• w.  -terminal  negative-resistance  devices  as  well  as  three-terminal 

.  A  symmetrical  fundamental  and  harmonic  power-combining  scheme 


I 


I 


was  developed  that  should  be  useful  in  extending  the  frequency  of 
operation  of  solid-state  devices. 

>.  The  important  elements,  such  as  FETs  and  varactors  that  arc- 
needed  for  implementation  of  the  symmetric  power-combining  scheme  in 
a  monolithic  integrated  circuit  form,  were  fabricated  and  characterized. 

6.  A  theoretical  analysis  was  carried  out  on  a  finite  YIG  slab 
in  a  rectangular  waveguide  and  for  properties  of  the  magnetostatic- 
wave  propagation  inodes  and  their  dependence  on  the  YIG  slab  dimensions 
and  its  location  as  well  as  the  magnetic  field  were  determined. 
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